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Abstract 

In a recent paper, one of us studied spherically symmetric, asymptotically flat solutions 
of Shape Dynamics, finding that the spatial metric has characteristics of a wormhole - 
two asymptotically flat ends and a minimal-area sphere, or ‘throat’, in between. In this 
paper we investigate whether that solution can emerge as a result of gravitational collapse 
of matter. With this goal, we study the simplest kind of spherically-symmetric matter: an 
infinitely-thin shell of dust. Our system can be understood as a model of a star accreting a 
thin layer of matter. We solve the dynamics of the shell exactly and find that, indeed, as 
it collapses, the shell leaves in its wake the wormhole metric. In the maximal-slicing time 
we use for asymptotically flat solutions, the shell only approaches the throat asymptotically 
and does not cross it in a finite amount of time (as measured by a clock ‘at infinity’). This 
leaves open the possibility that a more realistic cosmological solution of Shape Dynamics 
might see this crossing happening in a finite amount of time (as measured by the change of 
relational/shape degrees of freedom). 


1 Introduction 

Shape Dynamics (SD) is a Hamiltonian theory which describes gravity as the evolution of a 
3D conformal geometry. In this description, 4D spacetime is not the fundamental dynamical 
object of the theory, and must be understood as an emergent concept. Nonetheless it can be 
useful to describe the behaviour of weakly-backreacting ‘probes’. The utility of the space-time 
picture comes about from the relation of the two theories: SD is equivalent to the Hamiltonian 
formulation of General Relativity (GR) [1] for solutions of the latter which possess a complete 
‘constant mean curvature’ foliation. This slicing condition fails at the global level, already for 
simple examples: Schwarzschild’s spacetime does not possess a complete CMG slicing. In [2], one 
of us found that, respecting spherical symmetry and standard asymptotic flatness conditions, a 
solution of the equations of SD could be built which covers only the two non-singular quadrants 
of the Kruskal extension of Schwarzschild. The spatial metric is that of a wormhole, possessing 
a ‘throat ’(a minimal-area sphere) and two asymptotic ends. At the throat the spatial metric is 
smooth, but the spacetime metric has a defect. This discrepancy is possible in spite of Birkhoff’s 
theorem, because in SD it is the spatial conformal geometries that need to be regular, whereas 
in GR it is the 4D spacetime geometry that needs to be well-bahaved. 

The result of [2] was nonetheless preliminary: the assumption of spherical symmetry leaves 
no room for degrees of freedom of the conformal geometry (the Gotten tensor vanishes [3]). 
Moreover [2] did not address the issue of whether this solution would form from gravitational 
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collapse, i.e. whether it was physical or not. Here we partially remedy both shortcomings in 
a minimal way, by considering spherically symmetric matter sources, which introduce genuine 
dynamical degrees of freedom for the shell, ^ and studying their gravitational collapse. 

As was done in most recent works in SD [2, 3, 4], we exploit the local equivalence between 
SD and GR in CMC foliation to simplify the calculations. The input from SD is limited to; 
i) insisting that the spatial conformal geometries remain regular throughout evolution, and ii) 
neglecting any regularity requirement on the spacetime metric. 

There are other flaws in [2] which we are also not going to address in this paper. The first 
is the assumption of asymptotic flatness; according to the relational underpinnings of SD, the 
theory should be based on a closed spatial manifold. The second is the boundary conditions 
at infinity, which have been borrowed wholesale from GR, while in SD they should be set by 
the physical behaviour of matter ‘at infinity’ [3]. The last one is the use of maximal-slicing 
time, which violates one of the relational pillars on which Shape Dynamics rests: that time 
should be derived from the change of physical (i.e. shape, or conformally-invariant) degrees of 
freedom [5]. Nonetheless, there are arguments that show that, as far as the solution exists, it 
will reliably represent a “background experienced” space-time for weakly back-reacting matter 
degrees of freedom [6], and maximal-slicing time should approximate the amount of change 
experienced by a clock far away from the origin. In sum, the goal of the present paper is to 
study whether the solution found in [2] can emerge as the result of gravitational collapse, and 
therefore we postpone addressing these issues to further studies. 

Differently from previous work, here we need to couple matter - pressureless dust in par¬ 
ticular - in a first attempt to model gravitational collapse. The simplest distribution of dust 
which respects spherical symmetry is an infinitely thin sphere. The coupling of SD to matter 
is borrowed (for phenomenological reasons) from GR, by working, again, in the gauge in which 
the two theories are equivalent. This is done in Section 3. Before that, we will solve the SD 
equations in vacuum in Section 2, and then insert those constructions where appropriate when 
treating the full system coupled to dust. The full system is described in Sec. 4, where the 
reduced phase space of physical degrees of freedom is characterized, and the on-shell orbits 
describing the evolution of the collapsing shell are found. Sec. 5 contains an outlook of the 
result. 

2 Vacuum spherically symmetric solutions 

The constraints of Shape Dynamics, in the gauge in which it is equivalent to maximal-slicing 
GR, are [7, 8, 5] 


^ {p'^Pij - - ^/9R, 'Ui =-2 Vjp>i, C=p, (1) 

where gij is the spatial metric, its conjugate momentum and p = gijp^^ ■ These constraints 
need to be valid at all times. The time evolution of the fields is generated by the following 
equations: 

9ij {Pij ~ 29ijP) + > 

(Rd _ yjR) + N ij _ 1^2) 

^ (2) 
- ^ [P^k^ - \ppy + Va (VV^iV - g^^AN) 

+^kip'^e) - f^^k^ - p'^^^kC , 

which depend on a Lapse function N and a shift vector The latter is fixed by our choice of 
coordinates, while the former is fixed by the requirement of preservation in time of the conformal 

^Albeit not genuine “shape” degrees of freedom: for those one wonld need to introduce more than one shell, 
to form scale-invariant ratios of their variables. 
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constraint C ~ 0. This gives rise to the so-called lapse fixing equation: 

\^NR-2^^N + = Q. (3) 

V9 

Using the notation of [3], the most generic spherically symmetric metric and momentum can 
be written as: 

gij = diag {g^,cr,a sin^ 0} = diag §, ^^^|sin 6 », (4) 

where g, a, f and s are functions of the radial coordinate r only. There is an analogue ansatz 
for the shift vector: = (^(r), 0 , 0 ). 


2.1 Solution of the constraints 

Replacing the ansatz (4) into the constraints (1), we get 

-|- 4/^/i^ — 4fag?s -|- 12aga" 

- Uaa'n' - - I2ag? = 0 (5) 

nf' — \sa' = 0 , /i/ -|- s(T = 0 , 

where ' refers to the radial derivative The last equation can be solved algebraically, s = —^ f, 
and after replacing this expression for s in (5), it is easy to see that the diffeomorphism constraint 
can be written as a total derivative, ifV^) = 0. The solution of this equation is 

/ = 4 . ( 6 ) 


where A is an integration constant (meaning that it is spatially constant but can, in principle, 
still be a function of time). Finally, with a little work one can check that the Hamiltonian 
constraint can be rewritten as 


o-g/j d 
g' dr 



4^4 



(7) 


where the term yf' + 5 / ^ j is identical to the diffeomorphism constraint and therefore vanishes 
on-shell. The remaining term is a total derivative, and we can solve the equation by introducing 
a new integration constant m (the —4 factor is introduced for convenience), 

> f2 

^ — 4^ - — = —8m. 



Replacing the solution ( 6 ), we get a relation between a and /i, and since the latter appears 
without derivatives, the easiest thing is to solve with respect to //: 


^2 


^ - 8my/a + 4a 


(9) 


We found a solution to all our constraints which apparently holds for any choice of the 
remaining free function (T{r). This is a reflection of radial diffeomorphism invariance. In fact 
notice how grr = is homogeneous of degree two in a': the expression g?(4r^ oc (u'dr)^ 
appearing in the metric is explicitly invariant under changes of radial coordinate. However the 
choice of cr(r) is not completely arbitrary. If we require regularity of the conformal geometry, 
there are obstructions to the values that a can take. In fact, by inspecting (9) we can see 
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how the right-hand side is not guaranteed to be positive. It relies on the following fourth-order 
polynomial of y = y/ajm 

?{x) = - 2 sign(m) , C = . (10) 

2 

See [3] for a detailed discussion of the roots of IP. Here we only need to observe that if (7 = 0 
T is positive when x > 2. 


2.2 Solution of the equations of motion 

The equations of motion require previous calculation of the lapse from Eq. (3). Under the 
assumption of spherical symmetry (which for a scalar function like the lapse is just N = N{r)), 
the lapse-fixing equation reduces to 


_4^_ 

\ /ifj cr^ 
/8u' 8a' 

— ( ~ H 

\ /i'^ fj,^a 


ifi'a' 4cj" 
+ 


[a' 


r\2 




fi'^a 




.r, 8N" 

N' + = 0 . 



( 11 ) 


The solution is then N 
solutions, which are 


Cl A^i + C 2 N 2 , a linear combination of the two linearly independent 

a' f 


Ni = 


' 2^V^J (a')2 


dr. 


( 12 ) 


where the symbol refers to the principal value integral, which is needed because its argument 
contains the term 




\a 


{a') 


r\2 


^ -8m^ + 4a] 


\ 3/2 ’ 


(13) 


which diverges when a approaches a zero of CP (which has to be an extremum of a [3]). This 
divergence has opposite sign on the two sides of the extremum (the left- and right- limite are 
opposite), and the degree of divergence is the same, so that the following quantity is finite: 


(cj')2 {a'y^Jf+^{a'y) 


(14) 


where r G (ri,r 2 ) is the point where a has its extremum. 

Once we have the lapse we can calculate the equations of motion for the metric, the first 
of Eqs.(2). Using the spherical symmetry ansatz we get that the geo and equations are 
identical, and completely fix the shift vector: 


^i = 6'■^ ifN + a)/a'. 


(15) 


Replacing the above solution of in the grr equation (as well as the solutions of the ADM 
constraints), we find that they depend nontrivially on the lapse. Eortunately, replacing the 
solution of the lapse fixing equation (12), we find that the two principal-value integrals never 
appear explicitly - they always have an r-derivative acting on them, and we can use the fact 


that 







[a‘ 


r\2 ’ 


d, 


/ 


^3/2^3d^ 


(a'r 


(16) 


to simplify the expression. The equation then reduces to 


— 96 m -p QA (^A + C2^ = 0 . 


(17) 
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In order for the above equation to hold for any choice of a{r) the only possibility is that 

C2 = —2A, 171 = 0. (18) 

We fixed one of the two integration constant present in the lapse, 02 - ci cannot be fixed because 
the system is reparametrization-invariant, which implies the freedom to specify the value of the 
lapse at a point. Moreover, we discovered a conserved quantity: the integration constant m. 
This quantity has the significance of ‘mass-energy’, and it is conserved because the system is 
spherically symmetric, so it cannot radiate its energy away in the form of gravitational waves, 
and a form of conservation of energy similar to what holds in field theory on Minkowsi spacetime 
holds. The quantity m is actually what the Misner-Sharp mass reduces to in vacuum [3, 9]. 

The equations of motion for the momenta, the second of Eq. (2), are identically satisfied if 
one imposes the conditions (18), and therefore add no further information. We have then been 
able to solve exactly both the constraint equations and the equations of motion in vacuum. This 
is a significant result that can be exploited to build dynamically meaningful solutions of Shape 
Dynamics, for example by using some localized spherically symmetric distribution of matter, 
which leaves most of space empty. 


3 Coupling a thin shell of dust 

A thin spherically-symmetric shell of dust is just the continuum limit of a homogeneous dis¬ 
tribution of particles on a sphere, each one of which move radially with the same speed. One 
can deduce the appropriate contribution of such a source to the constraints of Shape Dynamics 
from that of a point particle. 


3.1 New form of the constraints 


The Hamiltonian and diffeomorphism constraints of ADM gravity coupled to a massive point 
particle are 

n = <5(3) (x* - y*) V5*^' Pi Pj + ruo , 

T-Li = -y")Pi, 

where y® are the coordinates of the particle, pi its momentum and mo its rest mass. The above 
constraints can be straightforwardly derived from the standard Einstein-Hilbert action coupled 
to a point particle. Note that pi is included as a cotangent vector, and this arises from minimal 
coupling. It is not hard to show that the constraints above are first-class. 

Now take a uniform distribution of point particles on the surface of a sphere of radius Rs, 
and take the continuum limit. The constraints become 


n = Vhp{Rs) 5{r - Rs)\Jg'^^p^ -h mg , 
Ri = d'"iVhp{Rs) S{r - Rs)pr , 


where hab is the metric induced on the sphere by pij, and p{Rs) is a scalar function to be 
determined [without weight: the additional weight 1/3 is provided by the delta function 5{r — 

Rs)]- 

To determine p{Rs) we have to ask that changing the radius of the sphere does not change 
the number of particles n: 


J dOdcpdrVh p6{r — Rs) = p{Rs) 


47rn, (21) 


which fixes p{Rs)- Now we can rescale the momentum npr = Ps, and the rest mass nrriQ = M 
of the single particle into the momentum and the mass of the whole shell, so that n drops out 


5 





of the equations. Now we can integrate over dOdcf), 


j nd9d(j) =-4Tr5{r - Rs)y^^rP[T^, 
j Tii dddcj) = Att S^i6{r — Rg) Pg , 

and the three constraints (5), are modified by the addition of the thin shell into 
- + 4/2//^ - AfaiJ?s + 12cr/rcr" - 12aa'fi' 

- 3^l{a'f - 12a^l^] = 6{r - Rs)^^+M\ 

/i/' - ^sa' = -^S{r - Rg), /r/ + su = 0. 


( 22 ) 


(23) 


3.2 Jump conditions 

After solving wrt s the maximal-slicing constraint, s = —^f, we can rewrite the second con¬ 
straint in (23) as 

= -lPsS{r-Rg). (24) 

The above equation has the form 

F'(r) = G{r) 6{r - ro) = G(ro) d{r - ro), (25) 

in any open set which does not include vq the solution to such an equation is F{r) = const. But 
because of the delta function on the right-hand side we cannot assume the continuity of F. In 
fact, integrating the equation from r = 0 to r one gets 


where 0(x) 


F{r) — F{0) = G{ro) 0(r — tq) -|- const ., 

0 , X < 0 

is the Heaviside distribution. Alternatively we can write 

1, X > 0 


(26) 


F{r) = F_ 0(ro — r) -|- F+ 0(r — ro), 
F+ - F_ = Giro). 


So (24) is solved by 


f^r) = ^eiRg - r) + ^0(r - Rg) 


a 


a 


with the ‘jump condition’ 


A I - A_ = 


2/r(i?s) 


(28) 


(29) 


For the Hamiltonian constraint, we should look at Eq. (23) and check which terms on the 
left-hand side can be divergent at r = i?s- / and s are not derived, and therefore they can 
contribute at most with a theta function, but not give any Dirac delta, /x and a have to be 
continuous because they are components of the metric. Therefore their first derivatives and 
a' can be at most discontinuous but not divergent, like / and s. The second derivatives 
and a", however, can be divergent, if the first derivatives are discontinuous. The only second 
derivative that appears is that of cr, so we can write 


singular part = —(5(r — Rs)^J 


(30) 
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So we have to assume that a is continuous, its first derivative has a jump, and its second 
derivative produces a Dirac delta term. The prototype of one such function is 


y{r) = yi(r) + {y^iv) - y 2 {Rs))&{r - Rs ), 
where yi{r) and y 2 {'>’) are continuous functions. Taking its second derivative: 

y” =y'[ + 2/20(^ - ^s) + 2?/2<J(?’ - Rs) 

+ (y2 - y2{Rs))S\r - Rs) , 


(31) 


(32) 


a distribution of the form z{r)6'{r — Rs) is not simply z{Rs)S'{r — Rs): with a smearing it is 
easy to show that it is equivalent to the distribution —z'(r)d(r — Rs) + z{Rs)5'{r — Rs), and in 
our case z{Rs) = limj._s.R 3 ( 2 / 2 (?”) “ V^iRs)) = 0. Then the above equation reads 


y" = y'i + y 2 ir)Q{r - Rs) + y' 2 iRs)S{r - Rs). 


(33) 


In terms of y, it is easy to see that the divergent term in y”{r) can be written as y 2 {Rs)S{r—Rs) = 
^lim^^^+ y'{r) — lim^_^^- y'{r)'^ 6{r — Rs). Then the jump condition for a" can be written as 


a” = lim (j'{r) — lim cr^(r) ] (5(r — i?s) + regular part. 




r^Ra 


(34) 


This (times 2/p) is the only divergent part of Eq. (30), and therefore we can identify it with 
the right-hand side: 

lim cr'— lim a' = + M‘^p?{Rs ), (35) 


r^Ra 


1 — >Rc. 


to produce our second jump condition. 

It is convenient, at this point, to define some quantities which will appear in all the jump 
conditions below: 

lim a = j, lim a' = k , cr{Rs) = , (36) 


r^R+ 


1 — >Rc. 


whatever diffeo gauge we choose, around r = Rs, a{r) can be written as 

a =p^ + {r - i?s )7 - Rs) 

+ {r- Rs)kQ{Rs -r) + 0[{r- i?s)'] , 

then, from the expression above, it is easy also to deduce that 

lim a = 2pp — 'yRs , lim a = 2pp — kRs . 


r^R+ 


r^R„ 


(37) 


(38) 


We are now in position to demand the continuity of p. Its expressions inside and outside of 
the shell do not coincide: 


A2 


, r < Rs 

W\/\l r>Rs 

and we have to demand that the left and right limit of p coincide: 

lim p{r) = lim p{r), 


r^R+ 


r^R„ 


that is. 




I 7 I 


A_ 

2p2 




(39) 


(40) 


(41) 
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This is a new equation we have to take into account, together with the jump conditions above, 
which, in the new notation, can be written as 


^-K = -yPi + Mif,^Rs). (42) 

We can completely eliminate 7 and k from the two jump conditions (41) and (42): using 
Eq. (29) into Eq. (42) 


K 7 

\^i{Rs)\ ~ \KRs)\ 


V , 


(43) 


by taking twice the square of the above equation, we can make it independent of the signs of k 
and 7 , 


+ 


7 


(^+- 41 -) 


Now, using the definition of fi{r) at r = Rs'. 


pI 



p\Rs) ■ 


(44) 


^ TO) = + V 

we end up with the following ‘on-shell condition’; 

- 4 ( m ++ m _)/9 + 4/)2 - = 

— 8 p + Ap^^ — 8 m- p + Ap^^ . 


(45) 


(46) 


3.3 Symplectic structure 

In order to discuss the dynamics of the system, we need to know which of the reduced-phase- 
space variables are canonically conjugate to each other. In other words, we need to calculate 
the symplectic form. By definition, the conjugate variables of the extended phase space are gij 
and , as well as Rs and Ps ■ Therefore the pre-symplectic potential is 


6 = J dr dO d(pp^^ 6gij -|- Att PsSRs , 

restricting it through spherical symmetry and integrating in dOdcp: 

poo 

6 = An dr {2f 6p + s 6a) + An Ps6Rs 
Jo 


(47) 


(48) 


Now we may impose the maximal-slicing constraint pf = —sa, and the solution to the diffeo 
constraint (28), 


(49) 


9 = An J dr ^ 2 / 6p — ^^6aj + An PsSR^ 

=—An f dr—^6{f\Ja) + AnPs6Rs. 

Jo 
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now, using Eq. (29) we observe that the first term in the last equation cancels the second term: 


dvr [ dr^6{f^/a) 


70 

poo 

- Stt / -^6 [^_ 0(i?s -r) + A+ 0(r - Rs)] 

Jo v<^ 


= — Stt 


+ Stt 


M_ dr-^ + 5A+ r dr-^ 
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a 


'Rs 


{A+-A.)-$Si^6Rs 


'^cr{Rs, 

= - Stt dr^ + M+ H dr-^ 

Jo V Jrs V ^ 

-AnPsdRs. 


And therefore the symplectic potential reduces to 

6 = —Stt 


Jo V Jrs V . 


(50) 


(51) 


4 Thin shell in an asymptotically flat nonexpanding region 

So far we have kept everything as general as possible. We now need to specialize to a particular 
model by fixing some of the integration constants. 


4.1 Boundary conditions 

The manifold we are studying has two boundaries, one at r —>■ oo (asymptotic infinity) and one 
at the origin r —)• 0. At infinity, as we discussed in the introduction, we have to set A+ = 0, 
in order to be consistent with [2]. This might not be entirely physically justified from the 
perspective of SD [3], but but here we are focused on determining whether the gravitational 
collapse of our thin shell of dust can generate the solution of [2], and therefore we have to 
impose the same conditions at infinity. Moreover, it should be noted that these are the standard 
asymptotically-flat conditions in GR [10, 11]. 

Regarding the inside of the shell, we cannot assume A_ = 0 because that would trivialize the 
dynamics (if 4+ = 4_ = 0 then Pg = 0). But we cannot assume 4_ ^ 0 for the entire interior 
either, because the metric would then develop a singularity, or a ‘piercing’-like defect (see [3]). 
We have then to assume that there is some other matter inside the shell, whose expansion 
compensates 4_ and puts the effective value of the integration constant A at the origin to 
zero. A realistic model of such matter could be, for example, a homogeneous-density star which 
accretes our thin shell as an additional layer. The price to pay is that m_ inside the shell cannot 
be put to zero either, otherwise we would be assuming the existence of some kind of matter 
with nonzero momentum but vanishing mass-energy. We thus set m_ as a free parameter of our 
model. We will ignore the dynamics of this conjectured matter near the origin, and concentrate 
on the exterior of the shell. The integration constants m+ and m_ are conserved quantities 
which characterize our solutions as freely adjustable parameters of the model. The only thing 
we assume about them is their positivity, because they are related to the Misner-Sharp mass 
and its positivity follows from the dominant energy condition when 4+ = 0 [3].^ 

^Although the condition is borrowed from the space-time picture, it should be valid in the appropriate limits. 


9 














4.2 Phase space 

Eq. (51), using the isotropic gauge /i = y/o jr, becomes 

9 = Stt log Rs S (4._|_ — A-) + exact form , (52) 

and, recalling Eq. (29), A^—A- = —^Ps Rs, we get (modulo an exact form) 9 = —AttPs logics ^Rs— 
AttRs logiis ^Ps, which gives the following canonical symplectic form: 

uj = 69 = A7r 6Ps A 6Rs , (53) 

so, in this gauge, Ps and Rs are canonically conjugate. 


4.3 Exact solution of the constraints 


So we set 4.+ = 0. Then in the region outside the shell we can use ‘isotropic’ coordinates, 
= -^, so that the metric outside is conformal to the Euclidean metric: 


a 

^ ’ 


ds^ = [dr^ + [d9^ + sin^ 6*d(/>^)] . 

Using this gauge, Eq. (9) can be treated as a differential equation for a: 

{ar 

j^2 

— 8 + 4(7 

and, using 4_|_ = 0, we can integrate this equation as 

/ 2,/a — 2 m_|_ + 2,/a — 2m^,/a\ 

I 2 m+ A:+ I 

where is a positive integration constant. Solving for a: 


2 m+ 


±1 


mz 


a = 


2 m+ 


±1 


+ 


1 




2 m+ 


±1 


(54) 


(55) 


(56) 


(57) 


an explicit calculation immediately shows that the above expression is identical whichever sign 
we choose (modulo a transformation —>■ 1/k^), so we can write 


mt 


a = 


A:+ r 
2 m+ 


+ 


1-1 4 

2 m+^ 2 
k+ r 


(58) 


The minimum of a is always 4m^, which is where \/o'/m+ reaches the only zero of the poly¬ 
nomial 7{y/a This minimum is at the coordinate radius r = 2m+/A:+. The integration 
constant has the only role of rescaling the coordinate r by a constant factor, and it is there¬ 
fore an effect of a residual radial diffeomorphism redundancy. We can fix this redundancy by 
imposing a - > r^, which means A:+ = 4. 

r^oo 


4.4 Solution of the jump conditions 

The solution (58) is valid outside of the shell. Inside the shell a will be different, because A- / 0. 
However in this region we cannot analytically solve Eq. (9) in isotropic gauge. Whatever the 
solution turns out to be, it will depend on one integration constant /c_. We need to satisfy two 
conditions, Eq. (41), imposing the continuity of /r (which, because we are working in isotropic 
gauge, implies also the continuity of a), and Eq. (42). The two equations depend on the left- 
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and right- derivatives of a at the shell, k and 7 , which are in turn determined by the two 
integration constants k- and k+. We can assume that k- has been solved by one of the two 
equations (41) and (42), and the other independent condition will be Eq. (46), which, after 
imposing A+ = 0 and A- = Ps (valid in isotropic gauge), is: 

(—4(m+ -I- m_) p + — |M|)^ = 

{4p‘^ - 8m+p) - 8 m-p + Ap"^^ . ^ ^ 

The equation above depends on through p = sja{Rs). To further reduce the number of 
parameters, we express everything in units of m+: 


Rs = m+ i? , Ps = m+ P , m_ = m+ a , Ms = m+ M , 


then Eq. (59) becomes 


^ {2R+1A 

64 W 



1)2 - m2 


16(2fl-|-l)2 


(l-2R)^R^ „2 

(2R+iy^ ^ 


(60) 


(61) 


In Eig. 4.4 we plot the on-shell curves P vs. R, for any possible choice of rest-mass M, and for 
a set of choices of a. Notice that the constant a, on physical grounds, should be smaller than 
one (and larger than zero), as the ADM mass inside the shell should be smaller than outside. 

We conclude this section with an analysis of the Hamiltonian vector flow in reduced phase- 
space {Rs,Ps). Consider Eq. (59) as a condition on if we take into account the def¬ 

inition of p = ^/a{Rs), it turns into an eight-order equation for m_|_. Let us write it as 
E(m_|_, Ps; M, m_) = 0. Its solution gives the ADM energy m+ as a function of the dynami¬ 
cal variables Rs and Ps (as well as the constant parameters m_ and M). This is the Hamiltonian 
generator of evolution in maximal-slicing time. If all we are interested in are the equations of 
motion of Rs and Ps in this time variable, we can avoid having to explicitly solve P = 0. We 


can instead differentiate F wrt all of the dynamical variables: -^^dmj^ + ^^dRs + ^^dPs = 0, 


dF 


dF 


which implies that ^ = - 
equations of motion generated by m_|_ are 


_ 


and -r^ 

F=o 


— _ dF. / dF 
~ dPs ! dP+ 


F=0 


Then the Hamiltonian 


Rs 


OF ( dF\ ^ 
dPs \dP+) 


Ps = 


^ ]ny_ (jFFy '' ^ 


dRs y dmj^ 


F=0 


(62) 


which, before replacing the solution P = 0, are two perfectly tractable functions of Rs, Ps and 
m+. 

Eq. (62) allows us to study the Hamiltonian vector ffow in phase space. In particular we 
can check its behaviour at the ‘throat’ Rs —?• m^j^l without having to solve P = 0. It turns out 
that the vector ffow vanishes at the throat: 


Rs -^ 0 , Ps-> 0 . (63) 

i?S—>-m+/2 iJs^m_|_/2 

The limits before are the same irrespective of the direction they are taken from, so, as expected, 
in maximal-slicing time the shell ‘freezes’ at the throat. One can also prove that it takes an 
infinite amount of maximal-slicing time for the shell to reach the throat, by explicitly integrating 
the vector flow. Maximal-slicing time has no intrinsic physical meaning: one of the fundamental 
relational underpinnings of SD is that time should be abstracted from the change of physical 
(i.e. shape) degrees of freedom. In this sense maximal slicing time is associated to the change 
in the DOEs of a clock far away from the origin. 


II 



















a = 0 


a = 0.4 



Figure 1: Plots of the on-shell surface (61) in the space R G [0, oo), P G (— 00 , 00 ), M G [0,oo), 
and for certain fixed values of a = G [0,1). The three variables R, P and M have been 

compactified by taking their arctan. The red curves represent the constant-M cross-sections, 
which are on-shell curves in the phase space R, P. Notice how all the curves ‘bounce’ on the 
P = ±00 boundary of phase space at R = 1/2. 

5 Outlook and Conclusions 

We are now in position to give at least partial answers to the questions we set forth at the 
beginning. First, does the ‘wormhole’-like line element found in [2] emerge from the gravitational 
collapse of spherically symmetric matter? Under the same assumptions of asymptotically flat 
boundary conditions (i.e. —?• = 0) at infinity the answer is clearly positive. The 

line element given by the areal radius (58) outside of the shell when = 4 is identical to that 
of ref. [2], so, as it collapses, the shell leaves in its wake the wormhole line element. 

The ‘on-shell’ relation (61) produces, for any value of M = and a = //ff a curve in the P- 
R space, which reaches the boundary of phase space P —)■ ±00 at R= that is, Rs = This 
value of Rs coincides with the throat of the wormhole line element with mass m_|_. This result 
implies that the collapsing shell does not reach the throat in a finite maximal-slicing time. This 
time parameter coincides with the experienced reading of a clock of an inertial observer at infin¬ 
ity. The preliminary conclusion is that the shell ‘freezes’ at the throat and cannot be observed 
to cross it. However, as we know, maximal-slicing time can at best be an infinitely-thin layer of 
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York time (the time parameter of CMC slicings). Whether the shell crosses the throat or not thus 
has to be postponed for the study of gravitational collapse in a cosmological setting (in which we 
take into account a cosmological constant, a nonzero York time and a compact spatial manifold). 

If this behavior does arise as a limit of the cosmological setting, at this point we would offer 
a tentative interpretation: the ratios of scales in a closed space-time (total volume, cosmological 
constant, and MS mass) may only allow for a given (non-zero) minimum areal radius (as it only 
allows for a maximum one). Thus either the system undergoes a bounce, or the shape degrees 
of freedom around the throat asymptotically (in time) freeze with respect to other local shape 
degrees of freedom. In either case, the dynamical behavior seems to be non-singular. 

For the moment, we can study the on-shell curves of (61), and observe that they continue past 
the point R = ^ where they reach the boundary of phase space. The solution curves fall into two 
topologically-distinct kinds: the closed and the open ones. The former are closed loops which 
touch the boundary of phase space at two points. They correspond to the cases in which the 
shell does not have enough kinetic energy to reach infinity, and recollapses back. Interestingly, 
this behaviour is observed on both sides of the throat R= so the shell recollapses also when 
it is in the region beyond the throat. The other kind of curves are the open ones, which reach 
the asymptotic boundary i? —)■ oo, and the other asymptotic infinity at i? —)■ 0 . 
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